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We present a fast and quasi-deterministic protocol for the production of single ions and electrons 
from a cloud of laser cooled atoms. The approach is based on a two-step process where first a single 
Rydberg atom is photo-excited from a dipole-blockade configuration and subsequently ionized by 
an electric field pulse. We theoretically describe these excitation-ionization cycles via dynamical 
quantum maps and observe a rich behavior of the ionization dynamics as a function of laser Rabi 
frequency, pulse duration and particle number. Our results show that a fast sequential heralded 
production of single charges particles is achievable even from an unstructured and fluctuating atomic 
ensemble. 



PACS numbers: 81.15.Jj, 32.80.Rm, 42.50.Hz 

The ability to place single ions into a medium or onto 
a surface with high precision opens up exciting possibili- 
ties for new types of nanofabricated devices and processes 
in materials science [THl]. Techniques such as scanning 
tunneling microscopy (STM) [B] and focused ion beam 
single ion implantation (FIBSII) [7] have been very suc- 
cessful in accurate placing of single ions. However, STM 
is relatively slow and FIBSII is stochastic and therefore 
governed by Poisson statistics. Ideally one would like 
fast, precise and fully deterministic single ion delivery. A 
promising route towards this goal is to use laser cooling 
and trapping techniques [8J, which allow exquisite con- 
trol over neutral atoms, therefore enabling 'atom optics' 
to become a realistic prospect for nanotechnology appli- 
cations [9]. By photoionizing ultra-cold atoms one can 
then transfer the precision control of neutral atoms on 
to charged particle sources. This approach has sparked 
research into the development of a new generation of 
monochromatic ion and electron beams pI)l - fT2"] . 

Thus far there have been two main approaches to mak- 
ing single ion sources based upon laser cooling. Firstly, 
laser cooled ions held in a trap can be ejected on de- 
mand. For example, linear Paul traps can hold strings 
of ions which can be emitted deterministically [131 114] . 
Another approach is to trap a single atom, for example in 
a magneto-optical trap (MOT) [TS] or optical trap [IB] , 
and then photoionize it. Whilst these approaches show 
great promise, the trap loading remains random and at 
present the repetition rate is slow (of the order Hz). 

In this work we propose an approach that bypasses 
the relatively slow process of carefully preparing an ini- 
tial state. This enables the fast and quasi-deterministic 
sequential production of ions out of a standard sample 
of laser-cooled neutral atoms. By analyzing in detail the 
quantum dynamics of the ion source we identify various 
dynamical regimes and discuss the temporal ion emission 
statistics as a function of experimental parameters and 



atom diffusion. We show that the system is not only 
suited for technological applications but also provides an 
interesting platform for the study of non-equilibrium dy- 
namics. 

The protocol underlying the envisioned ion source is 
based on a two-step process: At first, a tightly focussed 
laser beam that is resonant with the transition to a high- 
lying Rydberg state is irradiated onto an ensemble of 
ground state atoms. The use of Rydberg states ensures 
that at most one ion is created in the subsequent step, 
which consists of an electric field pulse ionizing the Ry- 
dberg atom. The reason for the excitation of only a sin- 
gle atom is the dipole-blockade [PTHST] which is a conse- 
quence of strong and long-range dispersive forces acting 
among atoms in Rydberg states. These lead to energy 
level shifts that prevent the laser-excitation of many- 
body states containing more than one Rydberg excita- 
tion, c.f. Fig. [lja). An experimental setup realizing the 
above scheme is sketched in Fig.[TJb). Here, an ensemble 
of ultra-cold atoms is prepared between two electrodes. 
The lasers exciting Rydberg atoms are focused such that 
the overlap of the beams creates an excitation volume 
with a characteristic size that is smaller than the blockade 
radius Rb- In a steady state MOT this can for instance 
be achieved by using a three-step excitation scheme |22) . 
The application of a small electric field results in the field 
ionization of the excited atom and the resulting ion and 
its concomitant electron are ejected from the ensemble 
in opposite directions as shown in the figure. The lighter 
electron can then be detected to 'herald' the creation of 
the ion. 

In what follows we will describe a simple model that 
captures the essential dynamical processes of such ion 
source and which will permit a detailed study of the 
statistics of the ion emission. Fig. [Ijc) shows the pulse 
sequence with which the ion source is operated. An 
excitation-ionization cycle consists of a resonant laser 
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FIG. 1. a) Dipole blockade. Energy level scheme of a pair 
two atoms that are resonantly laser-excited from the ground 
state \g) to a Rydberg state \r) with Rabi frequency Q. Due 
to strong interactions between Rydberg atoms the pair state 
\r,r) experiences a distance dependent energy shift. Below 
a critical distance Rb - the blockade radius - this shift be- 
comes larger than the Rabi frequency. Here \r,r) becomes 
energetically inaccessible and only a single Rydberg atoms is 
excited, b) Schematic of the ion source. Ultracold atoms are 
excited within the shaded volume through appropriately fo- 
cussed laser beams. The extension of the excitation volume 
is smaller than the blockade radius, Rb- Atoms can diffuse 
in and out of the region with rate Tin and r out respectively. 
An electric field is applied ionizing the Rydberg atom and 
a single ion and electron are emitted in opposite directions. 
The field generating electrodes (carrying charges ±ve) are in- 
dicated at the right and left hand side of the atomic cloud, 
c) Pulse sequence. An excitation-ionization cycle k consists 
of two pulses: a laser pulse of duration r ox that excites at 
most one atom to a Rydberg state (upper curve) and a field 
pulse of duration Ti on that ionizes the Rydberg atom (lower 
curve). The coherent Rydberg excitation is described by the 
atom number n dependent time evolution operator U n [see eq. 
|2|] and the ionization dynamics is determined by the Kraus 
operators M\ and M2 given in Eq. Q. 



pulse of duration r cx followed by an electric field pulse 
of duration T lon . The Hamiltonian describing the exci- 
tation step is given by H = ^ (a)b + b^aj, where f2 de- 
notes the single-atom Rabi-frequency, a* (a) a bosonic 
operator that creates (annihilates) a particle in the elec- 
tronic ground state and (b) an operator that creates 
(annihilates) a Rydberg atom. In order to take the Ry- 
dberg blockade into account in a simple way we im- 
pose fermionic anti-commutator relations {b, b^} = 1 and 
{b,b} = {b^tf} = on the b operators, so that b 2 = 0. 
This ensures that the excitation of two or more Rydberg 
atoms is not possible. 

The dynamics of the incoherent ionization step is gov- 
erned by a Markovian Master equation for the density 



matrix p of the (unionized) atoms contained in the exci- 
tation volume. In order to explicitly resolve the number 
of ions we project p onto the subspace in which m > 
ions have already been created. If p( m ) denotes the pro- 
jected density matrix then p — ^ m p^ m%> and the master 
equation governing its dynamics is [53] 
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with ionization rate 7. 

Since coherent laser-excitation and incoherent ioniza- 
tion happen in sequence, we can efficiently describe the 
dynamics of a whole excitation-ionization cycle within 
the framework of dynamical quantum maps 24 . The 
time evolution during the fc-th excitation-ionization cycle 
is given as follows: Starting at time tk-i the laser pulse 
coherently excites an atom in the excitation volume, i.e., 
the density matrix pk-i undergoes a unitary evolution 
p' = Upk-iU^ with U(t cx ) — exp(— iHr ox ). Since the 
laser-excitation does not change the number of particles 
in the volume, U is block diagonal in particle number 
space and the coherent dynamics for an ensemble of n 
atoms takes place in the subspace spanned by the states 
|p)® n and S [|s)® (n_1) ® \r)] /y/n, where S is the sym- 
metrization operator. The corresponding time evolution 
operator for the n particle block is 
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showing explicitly the y/n enhancement of the atom- 
field coupling due to the blockade effect [3T]. The sub- 
sequent ionization pulse maps the density matrix p' to 
p k = M\p/M{ + M 2 p'Ml, where 



\/l - e-T T »» 




(3) 



are Kraus operators. Mi describes the dynamics induced 
by the anti-commutator term of the master equation ([!]) , 
while M 2 captures the quantum jump processes induced 
by the first term on the right hand side of eq. ([T]). 

The density matrix of the subspace, where exactly 
m ions have been produced after the k-th excitation- 
ionization cycle can now be obtained iteratively 
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2 P k -! ' = M 2 U N _ m+1 pf-> W N _ m+1 Ml Here, 
denotes the number of atoms in the focal volume at t = 
and 0(x) the Heaviside step function. 

Let us now apply this theoretical framework to analyze 
the dynamics of the ion source in different operational 
regimes: 
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FIG. 2. Quasi-deterministic ion production. Probability p k 
of having produced k ions after k cycles for initial particle 
number N = 500 and excitation time t cx = ir/(yNH) (the 
other parameters are: Rabi frequency Q — 4f ionization pulse 
length Ti on = 167 -1 ). Inset: Probability distribution pi of 
having produced m ions after k cycles. The colormap is chosen 
such that white (black) corresponds to (1). The data shows 
that a single ion is created with nearly unit probability after 
each excitation-ionization cycle. 



(i) Quasi- deterministic regime in a microtrap - Here we 
consider an atomic ensemble which is confined within a 
trap volume which is much smaller than the blockade ra- 
dius. We start initially from a fixed number of trapped 
atoms N and the successive application of the excitation- 
ionization cycle will lead eventually to a depletion of the 
trap. For our numerical example we consider an initial 
atom number of N = 500 - a situation which can be 
achieved in case of an atomic Bose-Einstein condensate 
held in a tight optical trap. The Rydberg excitation is 
carried out with a 7r-pulse, with respect to the collective 
Rabi frequency y/Nfl. In the regime of a relatively large 
initial number of atoms, i.e. N ^> 1, this collective Rabi 
frequency changes negligibly with the emission of each 
ion. Thus, each excitation cycle of length t gx = w/y/NQ 
will - with probability close to unity - excite a single 
Rydberg atom, and this excitation process will remain 
efficient even after many ions have been emitted. This 
is indeed observed in Fig. [2] which shows the probabil- 
ity of having ionized exactly k atoms after k excitation- 
ionization cycles. In addition, the inset shows the proba- 
bility distribution of having ionized m atoms after k 
cycles. Hence, for a fully deterministic single ion source 
this quantity is 1 for m = k and elsewhere. The data 
indeed reflects that a single ion is produced in each cy- 
cle until the depletion of the trap begins to alter the 
collective Rabi frequency significantly. In principle the 
Rabi frequency could be adapted dynamically in order 
to account for atom loss. However, we can see that even 
in the static case, many atoms can be rapidly extracted 
from the sample with high probability. Using typical ex- 
perimental parameters, i.e. an excitation time r cx ~ 10 
ns and an ionization rate 7 ~ 100 kHz an ion production 
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FIG. 3. Trapped states. Probability distribution pi of hav- 
ing produced m ions after k excitation-ionization cycles for a 
fully blockaded sample initially prepared with N = 50 atoms, 
(a) The Rydberg excitation is performed using a single-atom 
7r-pulse, i.e., r cx = tt/Q.. (b) The Rydberg excitation pulse 
length is r cx = 4.927r/S7. The values of the other parameters 
are 7 = 4f7 and r; on = AO," 1 . The colormap is chosen such 
that white (black) corresponds to (1). 



rate up to 100 kHz appears feasible. 
(ii) Trapped states in a microtrap - We are now inter- 
ested in the ion emission dynamics in the case in which 
Rydberg atoms are excited with a single-atom 7r-pulse, 
i.e. r ex = 7ril _1 . The evolution of the corresponding 
probability density is shown in Fig. |3ja). We find that 
initially the mean number of emitted ions m grows with 
the number of cycles k. However, in this example with 
50 initial atoms, the probability of emitting ions beyond 
the 50-th cycle is strongly suppressed. In fact the data 
shows that no more than 14 ions will be emitted in to- 
tal even if the number of cycles is increased. The reason 
for that is the emergence of trapped states which occur 
when y/N — mf2T ex /2 becomes an integer multiple of tt. 
Here the time evolution operator |2]) becomes the iden- 
tity (up to a phase factor of ±1) and consequently the 
excitation step does not generate Rydberg atoms. This 
effect is analogous to the occurrence of trapped states 
in micromasers [251 126] . A further interesting dynamical 
regime is reached when the Rydberg excitation pulse is 
close to a v 7r-pulse with v = 1, 3, 5, .... In this case quasi- 
trapped states can occur that have a striking effect on the 
ionization dynamics. The corresponding probability dis- 
tribution for an excitation pulse length of r cx = 4.927r/f2 
is depicted in Fig. |3jb). As in the previous example a 
trapped state is reached after a few initial cycles. How- 
ever, since y/nQ.T eyi /2 can never be a strict multiple of it 
the trapping is not perfect and we observe a slow leakage 
into other quasi-trapped states. Eventually, this leads to 
a multi-modal ion distribution that peaks at ion numbers 
in where y/N — m fir ex /2 is close to a multiple integer of 

7T. 

(Hi) Excitation from an extended atomic cloud - So far 
we have considered a fully blockaded dense atomic sam- 
ple in a microtrap with radius much smaller than i?b- 
Let us now discuss the situation usually encountered in 
extended systems like standard MOTs. Here the atom 
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FIG. 4. Ion creation from an extended gas cloud. Probability 
pi of having produced m ions after k excitation-ionization 
cycles (a) without particle diffusion and (b) for diffusion rate 
K — 10~ 3 7- The initial number of particles in the excitation 
volume is distributed following a Poisson distribution with 
mean n = 10.2. The data shown are averaged over this initial 
particle number distribution. The excitation pulse length is 
T cx = 7r/vSf2. The Rabi frequency is f2 = 57 and the ion- 
ization pulse length is Ti D n = 207" 1 . The colormap is chosen 
such that white (black) corresponds to (1). (c) Influence 
of the choice of the excitation pulse length on the emission 
characteristics of the ion source. Shown is the probability 
of having produced k ions after k excitation-ionization cycles 
for t cx = ir/y/wQ (circles) and t cx = iv/Q, i.e. a single-atom 
7r-pulse (squares). All other parameters are those of panel 
(b). 



density is relatively small and, moreover, atoms can dif- 
fuse in and out of the excitation volume as depicted in 
Fig. [ijb) . In order to account for the effects of diffusion 
we use a simple rate model: Atoms can enter the excita- 
tion volume at a constant rate r; n and leave it at a rate 
Tout — Kn i where n is the number of atoms in the volume 
and k is the diffusion rate. In an equilibrium gas cloud 
we have r; n = r out and hence T m — nn where n is the 
mean number of atoms in the excitation volume. This 
simple model yields a Poissonian atom distribution with 
mean n for finding n atoms in the excitation region. 

To describe the influence of diffusion on the ion source 
dynamics we exploit the fact that the excitation cycle 
(r ox ~ 10ns) is typically much shorter than the ioniza- 
tion step (rion ~ 10/is). Thus we assume that diffusion 
is only relevant during the latter. This is consistent with 
the frozen gas approximation usually employed in the 
description of the laser excitation of cold Rydberg gases. 
Moreover, we consider the limit of slow diffusion, i.e. we 
assume KTj on 3> l/(2n), which means that during the ion- 
ization cycle only one atom will enter/leave the excitation 
volume at a time. Effectively this amounts to solving the 
diffusion equation up to first order. Within these approx- 
imations we can again use the above-employed theoreti- 
cal framework of dynamical quantum maps and arrive at 
dynamical equations analogous to Q (see Appendix). 

Let us now investigate the excitation dynamics. We as- 
sume that the initial mean number of atoms in the excita- 
tion volume is n = 10.2. With a realistic blockade radius 



of 5 /im this would correspond to typical MOT densities 
of 10 10 atoms/cm 3 . In Fig. [4] (a) and (b) we show the 
corresponding probability distributions f° r the case 
without diffusion and for a diffusion rate of k = 10 7, 
respectively. The data are averaged over the distribution 
of the initial number of particles in the volume and the 
excitation pulse length is fixed at t cx = ir/y/^tt. This 
choice therefore corresponds to a 7r-pulse with respect 
to the collective Rabi frequency \mfi. Comparing Figs. 
|4|a) and (b) we see that diffusion leads to a slow drift of 
the ion number distribution to larger mean values with 
increasing number of excitation-ionization cycles. This 
drift stems from the fact that the excitation volume is 
quickly depleted due the ionization process which is much 
faster than the diffusion. At long times (large k) further 
ions can only be produced when particles diffuse into the 
empty excitation volume leading to a slow increase of the 
total number of emitted ions. 

More interestingly, however, for the first few cycles the 
ion source produces one ion per cycle with a high prob- 
ability and is therefore quasi-deterministic. Note that 
this is despite averaging over the initial atom number 
distribution. Therefore, even when operated in an ex- 
tended atomic cloud, the behavior of the ion source is 
very similar to that in a fully blockaded sample discussed 
in (i). This is further corroborated by Fig. |4|c) where 
we show the probability of having emitted k ions after 
k excitation-ionization cycles. For comparison we also 
show the probability for that case in which the Rydberg 
excitation is performed with a 7r-pulse w.r.t. the sin- 
gle atom Rabi frequency fl. The data shows that the 
performance of the ion source is indeed significantly im- 
proving when the excitation pulse is chosen as 7r-pulse 
w.r.t. the collective Rabi frequency y/Ktl. Note, that the 
performance of the ion-source can be enhanced further 
by resetting the ion source after the few cycles during 
which quasi-deterministic ion emission is taking place. 
This is done by refocusing the excitation laser to a differ- 
ent part of the gas cloud thereby avoiding the depletion 
of the excitation volume and the broadening of the ion 
distribution visible in Fig. |4|a,b). 

In conclusion, we have presented and analyzed theo- 
retically a fast and heralded single ion source based on 
a dipole-blockaded atomic ensemble. The source can be 
operated in a variety of dynamical regimes and permits 
the quasi-continuous and quasi-deterministic production 
of single ions even from a extended gas cloud. The statis- 
tics of the ion emission shows some interesting dynam- 
ical features such as the emergence of 'trapped' states, 
reminiscent of micromaser physics. We expect that the 
presented scheme can be of practical relevance for appli- 
cations in materials science. 
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APPENDIX 



As discussed in the article we assume that diffusion is only relevant during the ionization step. Since we have two 
processes that can change the number of atoms in the focal volume (ionization and diffusion), we need an additional 
index to label the density matrix, p^"« m \ since the number of atoms in the excitation volume is no longer solely 
determined by the number of ions produced. Here, n refers to the number of atoms in the volume, m = 0, 1, . . . to 
the number of ionization events that have already taken place and a, f3 = (1) labels the electronic ground state 
(Rydberg state) of an atom. The evolution equations during the ionization step read 

&Voo' m) = 7Pn +1 ' m_1) - «(S + «)Poo'" l) + K ( n + 1 )Poo +1 ' m) + Kupfa' 1 '^ (5a) 

a,pir° = -7Pfr m) (5b) 

a tp ( r ) = _^+7 p ( r ) (5c) 

(n.m) K + 7 (n,m) / c 

dtPw = y~ p w . ( 5d ) 

where n denotes the mean atom number in the excitation volume at equilibrium (i.e., before starting to operate the 
ion source), 7 the ionization rate and k the diffusion rate. 

In order to obtain the dynamics of the excitation-ionization cycles we now need six operators: U n , M\ and Mi are 
those that we have already introduced for the diffusionless situation [eqns. (2) and (3) in the article] and in addition 
M 3 , Mf and M£ ut where 

t \ ( \A ~ KT(n + n) \ 
M 3 (n,T) = ( o x ) ( 6 ) 

and 

M?(r)=(^ Q \ , Mr { n,T)=(^ nQ \. (7) 



The operator M 3 gives the time-evolution of the density matrix that is diagonal in the number of particles n within the 
excitation volume, while the M4S describe the quantum jumps due to particle diffusion. The M4 operators change the 
number of particles in the volume without changing the number of ions to, while Mi affects both quantum numbers. 
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In the following it will be convenient to separate these two processes. To do so, we define block matrices that capture 
the dynamics within different n-subspaces by introducing a new set of Kraus operators. The block diagonal operators 
are 
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In contrast the jump processes are described by Kraus operators that have their non-trivial blocks on the sub- and 
super-diagonal in particle number space 
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Since the processes coupling different n-subspaces do not generate coherences between n and n' the density matrix is 
block diagonal in n 
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The time evolution of p^ m ' can now be described equivalently to the diffusionless case 

: M i (r ion ) M 3 (r ion )U{t cx ) p ( °l 1 U + (r ex ) M 3 (r ion ) At + (r ion ) 
+ Mi n r ion ) W(r ox ) p^\ U+ (r ex ) [Mf] + r ion ) 
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+ M^tw) W(r ox ) pj>> W+(T ex ) [X° ut ] + r ion )] (1 - 4 m ) + 
+ M 2 (rio„) W(r ox ) p^-V U+(r ex ) M+ (r ion ) J 0(fc - m) 
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(10b) 



These equations are the generalization of eqn. (4) in the article for the case that particles can diffuse in and out of 
the focal volume. The diffusion process is captured by the Kraus operators M 3 , M™ and M% ut . The operator A4 3 
captures the damping of the electronic ground state population as well as the electronic coherences between ground 
and Rydberg state for a given particle number n in the focal volume due to diffusion. The two M4 operators describe 
the (incoherent) change in particle number due to diffusion, which is treated as quantum jump processes. 



